Mon. Not. R. Astron. Soc. 000, [THED (2010) Printed 26 January 2013 (MN HT E X style file v2.2) 



Simulations of the formation and evolution of isolated 
dwarf galaxies — II. Angular momentum as a second 
parameter 

J. Schroyen 1 *, S. De Rijcke 1 , S. Valcke 1 !, A. Cloet-Osselaer 1 , H. Dejonghe 1 

1 Sterrenkundig Observatorium, Ghent University, Krijgslaan 281, S9, 9000 Gent, Belgium 
Accepted 2011 May 14. Received 2011 April 13; in original form 2010 December 14 

ABSTRACT 

We show results based on a large suite of N-Body/SPH simulations of isolated, flat 
dwarf galaxies, both rotating and non-rotating. The main goal is to investigate possi- 
ble mechanisms to explain the observed dichotomy in radial stellar metallicity profiles 
of dwarf galaxies: dwarf irregulars (dlrr) and flat, rotating dwarf ellipticals (dE) gen- 
erally possess fiat metallicity profiles, while rounder and non-rotating dEs show strong 
negative metallicity gradients. 

These simulations show that flattening by rotation is key to reproducing the ob- 
served characteristics of flat dwarf galaxies, proving particularly efficient in erasing 
metallicity gradients. We propose a "centrifugal barrier mechanism" as an alternative 
to the previously suggested "fountain mechanism" for explaining the flat metallicity 
profiles of dlrrs and flat, rotating dEs. While only flattening the dark-matter halo has 
little influence, the addition of angular momentum slows down the infall of gas, so 
that star formation (SF) and the ensuing feedback are less centrally concentrated, oc- 
curring galaxy- wide. Additionally, this leads to more continuous SFHs by preventing 
large-scale oscillations in the SFR ("breathing"), and creates low density holes in the 
ISM, in agreement with observations of dlrrs. 

Our general conclusion is that rotation has a significant influence on the evolution 
and appearance of dwarf galaxies, and we suggest angular momentum as a second 
parameter (after galaxy mass as the dominant parameter) in dwarf galaxy evolution. 
Angular momentum differentiates between SF modes, making our fast rotating models 
qualitatively resemble dlrrs, which does not seem possible without rotation. 

Key words: galaxies: dwarf - galaxies: evolution - galaxies: formation - methods: 
numerical. 



1 INTRODUCTION 

Morphologically, dwarfs come in two broad classes. 
Early-type dwarfs, or dwarf elliptical galaxies 
([Ferguson fe Binggelil 1 1994 dE;), are "red and dead" 
in the sense that their stellar populations are predomi- 
nantly old and that they are usually not actively forming 
stars. They almost completely lack the raw material for star 
formation: gas. In a small fraction of dEs, low level central 
star formation continu es at a rate of less than one solar 
mass every 1000 years (|De Riicke et al.ll2003al : iLisker et al. 
l200d ). dEs with luminosities below My ~ —14 mag are 
usually called dwarf spheroidals, or dSphs. As a class, dEs 
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are slowly rotating objects, flattened by velocity anisotropy 
(Ge ha et al.l [2003). L ate- type dwarf s, or dwarf irregular 
galaxies (dlrr; see e.g. ISkillman|[2005l V are gas-rich and are 
actively forming stars at a rate of about one solar mass 
every 100 — 1 000 years. As a class, dlrrs are flattened 
by rotation dCote et al.1 l200Ch . Noticeably, dSphs/dEs 
are found predominantly in dense galactic environments 
while dlrrs are typically found in more sparsely populated 
environm ents. This is the s o- called morpholog y- density 
relation dBinggeli et al.1 Il987l : ICote et al.1 120091 ). In the 
Perseus cluster, all dwarfs, irrespectiv e of type, appear t o 
avoid the very dense cluster center (| Penny et al.l [2009h . 
All this suggests that the environment is at least to some 
extent responsible for many of the differences between dlrrs 
and dEs. 

Despite their quite different properties, the two types 
of dwarfs also share many properties. They populate 
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roughly the same mass, m etallicity, luminosity, flattening 
(jBinggeli fc Popesculfl995h and size regimes and they have, 
to a good approximation, exponentially declining surface- 
brightness profiles. Moreover, the "boundaries" between the 
dwarf classes are not clear- cut and transition- type objects 
with mixed properties exist ([Grebel et al.ll2003h . This body 
of data provides us with evidence for evolutionary links be- 
tween, or at least a "common an cestry" for, the diff erent 
types of dwarfs. As shown by e.g. iMaver et al.l (|2006h . the 
combined action of ram-pressure stripping and tidal stirring 
on a star-forming, rotating late-type dwarf entering the halo 
of a Milky Way- like massive galaxy can remove most of its 
gas and angular momentum, effectively transforming it into 
a quiescent, non-rotating early-type dwarf. 

Dwarf galaxies entering a dense environment that 
are affected by ram-pressure stripping but not (or much 
less so) by tidal stirring, would be expected to keep 
many of their late-type structural properties and one 
would expect to find dlrr/dE transition- type dwarfs. In- 
deed, quiescent dwarfs have been observed that are signif- 
icantly more flattened and faster rotating than the aver- 
age dE (|De Riicke et al.1 l2003bl: Ivan Zee etaD l2004h con- 



tain gas and dust j Conselice et al. 2003; Buvl e et all 120051 : 
|Pe Looze et aljfeoiol ) , and often host embedded stellar disks 
and spiral structures (IJerien e t al. 2000; Ba razza et al ] |2002l : 
iDe Riicke erall l2003b: Graham et al. 200j). 



1.1 Metallicity profiles 

In general, dlrrs display chemical homogeneity practi- 
cally throughout their en t ire stellar and gaseous bodie s 
(jTolstov. Hill fe Tosil 120091 : iKobulnickv fe Skillmanl Il99/t) 
The Small Magell anic Cloud (S MCilDufour fe Harlowll 19771: 
Page! et al.1 Il978ft . N GC 6822 (iHernand ez-Mart mez et al.1 



An often quoted means of erasing metallicity gradi- 
ents in flattened dwarf galaxies is t he so-called "foun- 



20091 ). and Sextans A ([Kaufer et alj 12004 ) are examples of 



dlrrs without a significant chemi cal or abun dance gradient 
i n their gas content ; the SMC (|Cionil l2009h and IC 1613 
([Bernard et alJl20Q7h also lack a stellar metallicity gradient. 
Thus, a flat radial metallicity profile seems to be a rather 
general characteristic of dlrrs. iKoleva et al.l f|2QQ9h present 
radial stellar metallicity profiles, derived from optical VLT 
spectra, of a sample of 16 dEs belonging to the Fornax clus- 
ter and to nearby groups of galaxies. They find that ten of 
those, predominantly round and non-rotating, show a strong 
negative metallicity gradient. The six most flattened and 
most strongly rotating galaxies in the sample, however, show 
no significant gradient: like the rotationally flattened dlrrs, 
they are chemically homogeneous. Previous studies have also 
predominantly found neg ative metallicity gradients for dEs 
in the Local Group (e.g. lHarbeck et alJl200ll: lAlardl l200lL 
Sag gitarius; - the DART project: iTolstov et alJl2004 S culp- 
tor; [Battaglia et al hood. Fornax; Battaglia et al. 2010, Sex- 
tans), around M8 1 dLianou, Grebel, &; Kochll2010h and in 



the Coma cluster (den Brok et al.ll201lh . 



The findings of 



HP 

Koleva et al.l ([2009) suggest that, while 



total mass is most likely the dominant factor (as is concluded 

I 1 1 — I I 1 

from the simulations o f : iValcke et al. 2008; ISawala et alj 

l20ld : iRevaz et all l200d : IStinson et ai1l2007h . angular mo- 
mentum is an important second parameter in the chemical 
evolution of dwarf galaxies: fast rotating dwarf galaxies show 
a tendency to be chemically much more homogeneous than 
dwarfs with slow or no rotation. 
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references therein). The idea behind this mechanism is that 
the supernova feedback of a centralized star- format ion event 
is capable of ejecting significant amounts of hot, enriched 
gas through a cavity or "chimney" along the galaxy's minor 
axis. Subsequently, part of this gas can rain back down on 
the galaxy's disk, as in a fountain, diluting any metallicity 
gradient that might be present. In round low-mass galax- 
ies, centrally concentrated supernova feedback is expected 
to "blow away" all the gas rather than to "blow out" only 
the enriched hot gas. If this fountain mechanism is correct, 
the absence or presence of a metallicity gradient is deter- 
mined by two parameters: a dwarf galaxy's mass and its 
flattening or geometry. 

Alternatively, due to the "centrifugal barrier" in a ro- 
tating galaxy, gas cannot readily flow to the center and 
build up a strong centrally concentrated star formation 
event. One would therefore expect that rotation will nat- 
urally lead to more spatially extended star formation and 
thus to more spatially homogeneous stellar populations. In 
a similar vein, angular momentum has been proposed in the 
literature as the fundamental parameter setting low angular 
momentum starbursting Blue Compact Dwarfs apart from 
the m ore continuously star - forming high angular momentum 
dlrrs (Ivan Zee et alJl200lh . 



1.2 Paper 

In this paper, we use a suite of new N-body/SPH simulations 
to investigate how flattening affects the star- format ion histo- 
ries and chemical evolution of t h e isola ted dwarf galaxy mod- 
els presented in IValcke et"~ai1 (|2008h . Other N-body/SPH 
simulations of similar star-forming, gas-rich dwarf galaxy 
model s, though not always isolated, have been performed 
by e.g . lPelupessv, van der Werf & Ickd (l2004):IStinson et al.1 



by e.g . lremp' 

(1200a 120071) : IRevaz et al.l <|20Q9h : iGovernato et al.1 <|2010h : 
ISawala et al l (|20ld . l201lh . We flatten our originally spheri- 



cally symmetric models in different ways by adapting their 
initial conditions, with and without adding rotation, and 
compare the results both with the spherically symmetric 
originals and with the available observations. Our main goal 
is to contrast the "fountain mechanism" with the "centrifu- 
gal barrier" hypothesis, and to see if it's possible to produce 
dwarf galaxies with flat metallicity profiles in isolation. 

In section [2] we give a brief description of the simulation 
code, followed by a description of the simulations themselves 
in section [3] We present an analysis of the simulations in 
section 2] discuss the results in section [5] and conclude in 
section \6\ 



2 CODES 

For this research, we rel ied mainly on t wo codes: the Nbody- 
SPH code Gadget-2 (ISpringell [2005h for the simulations, 
and our own home-made analysis tool Hyplot for analysing 
and visualising these simulations. These will both briefly be 
described below. 
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2.1 Astrophysical mechanisms 

The code we actually use for our simulatio ns is a modifie d 
version of the Nbody-SPH code Gadget-2 (|SpringeJl2Q"05h . 
The freely available version only incorporates gravity and 
hydrodynamics, so to prove useful for investigating dwarf 
galaxy formation and evolution a number of additions were 
made. These include star formation, feedback and radiative 
cooling. Re-ionization or a UV background is not included 
in our models. Below we give a brief overview of the im- 
plementations, mo r e det ailed in formation can be found in 
IValcke et al.1 <|2Q08l , l20ld ) and in lValckd fcplOh . 



r m sNii, u 



Mi Mssp- 



Mi(m)$(m)dm 



2.1.1 Star formation 

Stars are formed when three criteria are satisfied: 

p g ^ pc — 0.1 cm -3 
T < T c = 15000K 
V.v ^ 0. 



(1) 
(2) 
(3) 



So we have a density treshold, a temperature treshold and 
the requirement that the gas be converging. We do not ex- 
plicitly implement a Jeans criterion. Gas particles eligible 
for star form ation are turn ed into stars according to the 
Schmidt law (|Schmidtlll959h : 



dp s 
dt 



dpg 

dt 



Pg 



(4) 



where p s , p g and c* are, respectively, the density of stars and 
gas, and a dimensionless star formation efficiency factor. t g 
is taken to be the dynamical time for the gas 1 / yj AnGp g . 



2.1.2 Feedback 

Produced star particles are represented as "simple stellar 
populations" (SSP), applying the Salpeter IMF for the prob- 
ability that a star of mass m resides in the SSP: 



$(m)dm = Am 



-(l+x) 



dm, 



(5) 



with x — 1.35, A = 0.06, and the limits for stellar masses 
are m\ = 0.01 M and m u = 60 M . Feedback from a star 
particle is given through stellar winds (SW) and supernovae 
(SN, type II and la), and includes the return of both en- 
ergy and enriched gas to the ISM (thermal feedback). These 
are transferred to the surrounding gas particles according to 
the SPH smoothing kernel. The lower limits for SNIa and 
SNII respectively are 3M and 8M , upper limits 8 M and 
6OM , and progenitor lifetimes are 5.4 x 10 6 yr and 1.5 x 10 9 
yr. 

The energy feedback for both types of SN is taken to be 
10 51 erg, for SW this is 10 50 erg, which are all transferred 
to the ISM with an efficiency of 0.1. 

The returned mass fractions for SNIa and SNII are 
0.00502 and 0.112, and the metal y ields from these super- 
novae are taken from, respe ctively, Travagl ic^eTaD {[2004) 
(their b20_3d_768 model) and lTsuiimoto et al.1 (|l995h . From 
the last authors we also adopt NsNia/NsNii = 0.15 to set 
the fraction of stars in the relevant mass range that reside 
in binary systems and that can go SNIa. The yield Mi of 
element i by e.g. SNII is then calculated as: 



J™ u m$(m)dm 



(6) 



The formula for SNIa is similar but simpler, since there is 
no dependence on progenitor mass (Mi(m) becomes Mi). 



2.1.3 Cooling 

Metal dependant gas c ooling is implemented accor ding to 
the cooling curves from Sutherland & Dopital (Il993h . These 
curves are interpolated in metallicity and temperature to ob- 
tain the cooling strength for a gas particle. They are allowed 
to cool in this fashion to a minimum of 10 4 K, below which 
further cooling is only possible through adiabatic expansion. 



2.2 Analysis and visualisation 

For the analysis we used our own Hyplot package, 
which is freely available on SourceForg^E It is an analy- 
sis/visualisation tool specially suited for Nbody-SPH simu- 
lations (currently only specifically for Gadget-2 datafiles), 
written mainly in Python and C++. PyQT and Matplotlib 
are used for the GUI and the plotting, and it is also fully 
script able in Python. All analysis, plots and visualisations 
in this paper have been made using Hyplot. 



3 SIMULATIONS 

In this section we will describe the simulations themselves: 
the models used for the basic initial conditions, the addi- 
tional setup we need for our goals, and finally the grid of 
our production runs, together with a preliminary evaluation 
of those runs to have an idea of the simulated objects we 
have at our disposal. 



3.1 Initial conditions 

We base the initial conditions of our flattened dwarf galaxy 
simul ations on the s pheric ally symmetric dwarf galaxy mod- 
els of IValcke et al.l (|2008h . We describe these briefly below 
and in Table [1] We introduce flattening into the models by 
adding initial flattening and/or rotation. The precise way in 
which this is done is described below. In the end, we have 
a set of flattened dwarf galaxies, both rotating and non- 
rotating. This way, we can distinguish between the effects 
of the geometry (flattening only) and the kinematics (i.e. 
rotation) . 



3.1.1 Spherically symmetric dwarf galaxy models 

The basic spherically sym metric dwarf galaxy models come 
from lValcke et al.1 <|2008h . They consist of 

(i) a dar k matter halo with a cored Ku z 'min Kutuzov den- 
sity profile ([Deionghe & de Zeeuwl ll988), 

(ii) a spherically symmetric homogeneous gas cloud, set 
to a fixed initial temperature of 10 4 K, initial metallicity of 
10 -4 Z© and the gas particles initially at rest. 



1 htt p : / / sour cef or ge . net / pr o j ect s / hyplot / 
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Table 



Details of the basic spherical dwarf galaxy models (see 
Initial masses for the DM halo and gas are in 



IValcke et al.ll20Q. 

units of 1O 6 M0, radius in kpc. 



model 


MnM,i 




a 


C01 


206 


44 


0.439 


C03 


330 


70 


0.513 


C05 


660 


140 


0.646 


C07 


1238 


262 


0.797 


C09 


2476 


524 


1.004 



We start off with only gas and dark matter (DM). The gas 
cools and collapses into the DM gravitational potential well. 
Our modified version of Gadget- 2, as described in the pre- 
vious section, then allows the gas particles to produce star 
particles if the criteria are satisfied (of which the density cri- 
terion proves to be the most important). Feedback of energy 
and of newly synthesized elements from these star particles 
is accounted for. For specific inf ormation and detail s on these 
models, we refer to TablefTland lValcke et all (|2008l ). The ac- 
tual initial condition files are created by random (Poisson) 
sampling of these density profiles for both components. 



constant 

arctangens 

combined arctangens/constant 
combined linear /constant 




Figure 1. The star formation rate as a function of time for four 
different rotation profiles (see text for details). 



is then kept constant at the constant value of the fiducial 
models ("combined linear/constant"). 

• the rotation velocity rises as an arctangens with radius 
out to 0.5 kpc and is then kept constant at the constant value 
of the fiducial models ("combined arctangens/constant"). 

These simulations show no noticeable influence of the central 
velocity profile on the behaviour of the models, see Fig. [1] 

We will henceforth refer to the set of CR models as the 
"rotating models". 



3.1.2 Adding initial rotation 

Rotation is added to the gas only in the initial conditions. 
The DM halo, implemented as a live halo, simply provides 
a background gravitational potential and is not given any 
rotation. Every gas particle is given a tangential velocity 
according to the desired rotation profile. We align the z-axis 
with the galaxy's rotation axis. 

Our preferred rotation profile is a constant radial rota- 
tion profile (CR), which means that the net velocity given 
to each gas particle is independent of radius. We do not use 
a solid body rotation profile (SBR), in which the velocity de- 
pends linearly on the radial distance from the rotation axis. 
This is because preliminary tests showed that using SBR, a 
significant fraction of the gas content of the galaxy imme- 
diately became unbound and was lost, even at low rotation 
speeds. The CR profile we use can be argued to have prob- 
lems on small radii. The constant value of the velocity means 
that the angular velocity rises quickly when approaching the 
rotation axis. But here our preliminary tests show that in 
practice this is not a problem. Because the gas particles 
are placed randomly and we are working with floating point 
numbers, no particle in practice will ever lie exactly on the 
rotation axis. The gas particles in the central region initially 
have a very high angular velocity, but interactions between 
the gas particles quickly slow them down and convert the ex- 
cess velocity to heat, wich is then very quickly lost through 
the very efficient radiative cooling. Also, since the gas is ini- 
tially distributed homogeneously over a large sphere, there 
is no large gas mass in this central region. We did a num- 
ber of test simulations with a variety of alternative rotation 
profiles with more well-behaved central velocities: 

• the rotation velocity rises as an arctangens with radius, 
with the asymptotic velocity equal to the constant value of 
the fiducial models. 

• the rotation velocity rises linearly out to 0.5 kpc and 



3.1.3 Adding initial flattening 

As a flattening parameter we consistently use the axis ratio 
q = c/a in this paper, where c is the shortened z axis and a 
is the axis in the x — y plane (we adopt axially symmetric 
models). Here, we flatten both the gas and DM distributions. 
The parameters of the Kuz'min Kutuzov profile naturally 
allow for introducing a flat tening to the DM ha lo in the form 
of an axial ratio (Deionghe AdeZee^3[l988), so the DM 
part of the flattening is trivial. Flattening the homogeneous 
gassphere takes a little bit more care. We need to scale up the 
axis in the x — y plane appropriately while scaling down the 
z axis to ensure that the density of the gas cloud remains the 
same. This is important because we do not w ant to change 
any a spect of the dwarf galaxy models from IValcke et al.l 
(2008), other than the geometry, and the initial density is 
a very important one. Starting from a sphere with radius 
r, if we want to achieve a flattening q we produce it by 
calculating the a and c axes as follows: 



This will ensure that the volume of the ellipsoid is constant 
for any value of q, and therefore so is the density. 

In the remainder, we will refer to this set of models 
as the "flattened models" (although of course the rotating 
models also become flattened eventually). 

3.2 Production runs 

In Table[2]we show an overview of the run numbers and spec- 
ifications of our production runs. We simulate dwarf galaxies 
with a range of masses, flatt enings and ro t ation speeds using 
the dwarf galaxy models of I Valcke et al] fl2QQ8h from Table 
[T] and additional methods for setting up the initial condi- 
tions as described above. Both the gas and the dark matter 
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Table 2. Grid of the production runs, given with runnumbers 
and specifications of the three used parameters: mass (first col- 
umn, see Table [T]) . initial flattening (q, last column) and initial 
rotation speed (vi, 3 different rotation speeds) (see above). 
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0.1 


C07 
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1 
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0.5 
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277 


287 


0.1 


C09 


209 


219 


229 


1 




239 


249 


259 


0.5 




269 


279 


289 


0.1 



components were represented by 200000 particles, and the 
simulations were evolved in time during 11.7 gigayears, cor- 
responding to the time from z = 4.3 to the present. The 
seeds used for sampling the particles from the specified den- 
sity profiles (see section I3.1.ip are chosen at random for all 
production runs. 



3.3 Preliminary evaluation of simulations 

Table \3\ lists many different physical quantities for all of 
the simulated dwarf galaxies in our set from Table [2j These 
are the final values for these quantities, evaluated at the 
end of the simulation, except those explicitly indexed with 
'i', which are initial values. Broad-band colours are cal- 
culated (with bilinea r interpolation) using the models of 
IVazdekis et al ] (|l996h . who provide mass/luminosity values 
for SSPs according to metallicity and age. For those sim- 
ulations that form little or no stars, making an accurate 
evaluation of the physical parameters impossible, we simply 
enter a "— " in Table [3] 

To evaluate our methods for setting up the initial con- 
ditions, we discuss the C05 models below (see Table [2]). 



3.3.1 Variance 

We first make note of the inherent variance in our mod- 
els. To this end we have produced a set of 25 simulations 
of the basic spherical C05 model, with different samplings 
of the dark matter halo and the gas sphere. For each sim- 
ulation, different random seeds are used to construct the 
initial condition. As shown by the 15.9th/84.1th percentile 



0.008 




t [Gyr] 

Figure 2. Depiction of the variance inherent to our models, us- 
ing a set of 25 differently sampled initial conditions. Upper panel 
shows the evolution of the star formation rate, lower panel shows 
the evolution of the total stellar mass. The black line is the mean 
curve for our set. The green band shows the area between the 
15.9th and 84.1st percentile (which are linearly interpolated be- 
tween the closest ranks, and would correspond to the la interval 
if the underlying distribution was Gaussian), and the light grey 
band shows the area between the minimum and maximum value 
of our set. These percentiles and extrema are calculated in each 
time-bin (of which there are 100). The two dark grey lines show 
the evolution of two individual runs: the dashed and dotted line 
represent the runs which, respectively, produced the lowest and 
the highest total stellar mass at the end of the simulation. 

area and the total range of the SFHs of this set of simu- 
lations in Fig. [21 the variance is significant, allowing for a 
variety of SFHs. This is however not unexpected. Systems of 
this kind, with stochastic star formation and feedback, are 
inherently chaotic. Small differences are continuously ampli- 
fied and can, over time, lead to large deviations. But on the 
other hand, more importantly, the green band depicting the 
15.9th/84.1th percentile area (which would correspond to 
the la interval if the underlying distribution was Gaussian) 
shows quite clearly the generic behaviour of the models. So 
we keep in mind that our models can show a spread in their 
properties, but that they also exhibit a clear general be- 
haviour. 

3.3.2 Rotating models 

As described above, to obtain a rotating galaxy we add ini- 
tial angular momentum to the gas. We need to check if 
this actually results in a rotating stellar component of the 
galaxy. In the upper panel of Fig. [3] we present the rotation 
curve of the gas particles at different times. In the lower 
panel, we show the final stellar rotation curve (at 11.7 Gyr). 
These are binned profiles of tangential velocity versus dis- 
tance to the z axis, where the profile value in every bin is 
the average rotation velocity per particle in that bin. 

The rotation profile of the gas rises due to the gas falling 
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Table 3. Details of simulations. All physical quantities are evaluated at the end of the simulation (11.7 Gyr), except those indexed with 'i', 
which are evaluated at the beginning. Columns: (1) model number (see Table [T}, (2) simulation number, (3) initial flattening (gas/DM), 
(4) initial rotation speed of gas [km/s], (5) spin parameter of gas in IC, (6) final gas mass [1O 6 M0], (7) stellar mass [1O 6 M0], (8) half- 
light radius [kpc], (9) (10) B-band and V-band magnitude, (11) (12) fitted Sersic parameters of surface brightness profile, (13) central stellar 
velocity dispersion along line of sight (edge-on) [km/s], (14)(15) luminosity-weighted metallicity (B-band), (16) final flattening of the stellar 
component (averaged over last 3 Gyr), (17) final stellar peak rotation speed [km/s]. Omitted values were irrelevant due to low stellar mass. 
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into the potential well, and quickly evolves to a rather stable 
form, only perturbed temporarily by the turbulence caused 
by strong star formation events. This "steady state" is a con- 
sequence of the balance between cooling, which makes the 
rotating gas sink inwards, and supernova feedback, which 



heats and disperses gas. The stars that form from the gas fi- 
nally follow a rotation profile that rises out to one half-light 
radius and flattens off beyond that radius. This confirms 
that using CR initial conditions for the gas is adequate to 
achieve stable, rotating dwarf galaxies. The final rotation 
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Figure 3. Rotation curves of our showcase model (225), upper 
panel displays the gas at different times during the simulation, 
lower panel displays the stars at 11.7 Gyr (with adaptive binning). 



speed of the stars which is included in Table [3] is the peak 
value of this rotation curve. 



3.3.3 Flattened models 

In Fig.Uwe show three of our model galaxies (a non-rotating 
spherically symmetric model on the left, a rotating model 
and a flattened model) to see the resulting flattenings at 
the end of the simulations. The latter two both show a con- 
siderable and stable flattening (see the C05 model in Fig. 
[5} . Some trends become apparent when looking at the total 
mass range in Table [3] and Fig. [5] It appears that the stabil- 
ity of q for the stellar component significantly increases with 
rising mass in the rotating models. Only in the least massive 
models does q rise significantly with time. The more massive 
models all exhibit a stable flattening around q ~ 0.4 — 0.5, 
so our model dwarf galaxies are relatively "thick". In non- 
rotating models with an initially strongly flattened halo, the 
halo thickens and it turns out to be impossible to make stel- 
lar bodies more flattened than 0.6 — 0.8. As can be read 
from Table [3l combining initial rotation and initial flatten- 
ing helps somewhat to achieve stronger flattenings in the 
least massive models. 

We thus note from Table [3] and Fig. [5] that we are 
not able to make ex t remel y flat galaxies. In the work of 
iRovchowdhurv et al a collection of dlrrs from the 

FIGGS survey is investigated, and they find from the 
flattening distribution a mean axial ra tio (q) » 0.6 for the 



HI dis ks. Si milar values are obtained by Binggeli &; Popescul 
(Il995h. IStavelev-Smith. Davies. &; Kinmanl (|l99^ L 
Hunter fc Elmegreenl <2006h. ISung et al. 1 (|l998h and 
Sanchez- Janssen, Mendez-Abreu, & Aguerr i (|2010h 
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Figure 4. Edge-on and face-on views of the stellar distributions 
of a non-rotating spherical model (left, 205), a rotating model 
(middle, 225) and a flattened, non-rotating model (right, 265). 
All are slices of thickness 0.4 kpc, axes are in kpc, and color 
denotes projected stellar density. 
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Figure 5. Evolution of the flattening parameter q = c/a of 
the stellar component during the simulations for different galaxy 
masses. Upper panel: rotating models (all with Vi = 5km/s), 
lower panel: flattened models (all with qi = 0.1). Per panel all 
properties are identical, except for the mass. Only galaxies with 
an appreciable stellar mass are shown (see Table [3]). Dotted lines 
show the q of the DM component, only shown in the bottom panel 
because the DM consistently has q = 1 in the top panel. 



for the stellar content of dwarf galaxies. Also other 
simulations suggest that low mass galaxies are not 
born as thin discs, but as thick, puffy systems (e.g. 
iKaufmann, Wheeler, fe Bullockl 120071 : and to a lesser 
exten t the (more massive) models of iGovernato et all 
I201CL which are still not extremely flat). The reason for 
this is sought in the increasing importance of turbulent 
motions, plausibly caused by star formation and feed- 
back, with respect to rotational motion in low mass systems 



(Kaufmann, Wheeler, & Bullock 2007; Rovchowdhury et al. 
l20ld : ISanchez-Janssen. Mendez-Abreu. & Aguerril l2010h . 
Besides the moderate value of the flattening itself, we 
also qualitatively reproduce the trend with gala x y mas s 
from ISanchez-Janssen, Mendez-Abreu, fc Aguerril (|2010h . 
All our simulations are below their 'limiting mass' of 
M* « 2 x 10 9 Mq, and indeed for simulations with identical 
initial setup, the final stellar bodies thicken with decreasing 
mass. This can be seen in Fig. [6] where we mimic their 
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Figure 6. Axis ratios of all our models. Left panel: versus stellar 
mass (in Mq), right panel: versus B-band luminosity. The differ- 
ent initial rotation speeds are indicated with color (white: km/s, 
grey: 1 km/s, black: 5 km/s), the initial flattenings are indicated 
with symbol shapes (lozenge: 1, triangle: 0.5, square: 0.1). 



Fig. 1 (the leftmost and rightmost panels) 
massive, flat models (M* « 3.4 x 10 8 M fl , 



Our most 
g g 0-4) 



conne ct nicely to models DG1 and DG2 of iGovernato et all 
(|2Q10h , who are slightly more massive and slightly natter 
(M* « 4.8 x 10 8 M , q « 0.35). 



between different models of different sizes, the half-light ra- 
dius of each simulation is also indicated on the plots with 
a dashed vertical line. We note that the rotating models 
can be considered to have fiat profiles out to 1.5 times R e L, 
while the spherical models usually show a fall-off well before 
that. Noticeably, the mean [Fe/H] of the lower mass models 
appears to be significantly higher when rotating, while the 
opposite is true for the higher mass models. This will be 
discussed further on in section 14.4.41 

Fig. [8] shows the same quantities for some of the flat- 
tened models. Surprisingly, there appears to be no obvious 
trend between the flattening and the shape of the metallicity 
profile, with most galaxies showing strong negative metal- 
licity gradients. We show all models from Table [2] which 
received an initial flattening but no initial rotation, so for 
each mass we have 2 different degrees of initial flattening 
(q = 0.5 and q = 0.1). It is clear that the flattening gen- 
erally has no significant effect on the metallicity gradient, 
almost all simulations have a negative slope. Only in the 
most massive ones, or where the spherical model does not 
have a strong gradient to begin with, does the initial flat- 
tening appear to have some ability to somewhat flatten the 
metallicity profile. 



3.3.4 Galaxy mass and concentration 

For a rotating model the half-light radius (R e L), defined 
as the radius of the sphere containing half of the light) is 
considerably larger than that of a sperical model, as can be 
seen m Tableland Fig.g] The total stellar mass usually de- 
creases when adding significant rotation. Non-rotating flat- 
tened models on the other hand are generally not much 
larger than the spherical models, sometimes even smaller. 
The half-light radius decreases slightly in flattened galax- 
ies at lower masses, and increases slightly at higher masses 
(with respect to the spherical model). The total stellar mass 
decreases with increasing flattening. 

The rotating models are thus spatially more extended 
than their spherical progenitor and at the same time they 
generally are also less massive (in stellar mass) so they 
are considerably less centrally concentrated. The flattenend, 
non-rotating models are usually less spatially extended than 
the spherical models and also less massive, so they have sim- 
ilar central concentrations (see also Fig. [4}. 



4 ANALYSIS 

In this section we present a more extensive analysis of our 
production runs. 

4.1 Metallicity profiles 

Looking at Fig. [71 we see some interesting results concerning 
the metallicity profiles of the rotating galaxies. For a range 
of galaxy masses we compare the metallicity profiles of the 
spherical models with those of the fastest rotating models 
from Table [2] (with v% — 5 km/s) in Fig. [7] The metallicity 
profiles of the spherical models almost always show a clear, 
negative gradient, while the profiles of the rotating mod- 
els are always significantly flatter. For a proper comparison 



4.2 Star formation histories 

Next we turn our attention to the star formation, and for this 
we look at Fig. [9] where star formation histories (SFH) of 
different simulations are shown. Rotation also seems to have 
a significant influence here. We again compare the spherical 
models with the fastest rotating models from Table [2] for a 
range of galaxy masses in Fig. [9j where we show the evo- 
lution of the produced stellar mass expressed in solar mass 
per year (M©/yr). The total stellar mass of the galaxy is 
plotted in dashed lines alongside the SFHs. 

Non-rotating spherical models typically have "breath- 
ing" or "bursty" SFHs, with strong SF peaks a few Gyr 
long, separated by quiescent periods wh ere the star forma- 
tion rate (SFR) essentially goes to ze ro dValcke et al. 1 120081 : 
IStinson etl!ll2007l : lRevaz et al.ll2009h . The strength and du- 
ration of these peaks, as well as the intermittent pauses, 
depend mainly on galaxy mass. 

The models with rotation, however, are able to reduce 
this burstiness and make the SFH much more continuous. 
Periods of increased star formation still exist, alternated 
with lulls, but the SFR never drops down to zero. The effec- 
tiveness of reducing the SF peaks varies in our simulations, 
and depends primarily on the galaxy mass. In the least mas- 
sive models, which in the non-rotating spherical case show 
the most extreme bursty behaviour (one big initial burst 
almost completely shutting down further SF activity), the 
effect of adding rotation is most noticeable. The SFR now 
becomes virtually flat. The more massive models still show 
some SF fluctuations but not nearly as pronounced as in 
their spherically symmetric analogs. 

Flattening on the other hand does not have a large effect 
on the star formation history of the galaxies. When looking 
at Fig. nja we can see that flattening, unlike rotation, gen- 
erally does not induce major qualitative differences in the 
SFH. The SFHs are generally very much like the SFHs of 
the spherical models, still having large peaks separated by 
periods with zero star formation. 
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Figure 7. Metallicity profiles of our productions runs, see Table 
[2] Each frame compares, for a certain galaxy mass, the spheri- 
cal model (20a;) with the fastest rotating model (22a?). The R e L 
of each model is also indicated with a dashed line, and for the 
rotating model we also show 1.5 x R e L for indicative purposes. 
Adaptive binning was used to produce these profiles, the width 
of each bin being indicated by a horizontal grey bar. 



Figure 8. Metallicity profiles of our productions runs, see Table 
[2] Each frame compares, for a certain galaxy mass, the spherical 
model (20a;) with the 2 non-rotating flattened models (23a: and 
26a;). Further details of the plot are similar to those of Fig. [71 



4.3 Gas structure 

The structure of the gas of dwarf galaxies is another typ- 
ical characteristic that we will consider. The best obser- 
vations available in this respect come from the Magellanic 
Clouds; for instance iKim et all (|2005h present HI data for 
the LMC. Another very useful source of observational data 
about the HI gas content and structure of dwarf galaxies 
is the THINGS survey (The HI Ne arby Galaxy Survey; 
IWalter et all 120081 : IWeisz et all l2QQ9h . These studies show 
that the neutral hydrogen gas of dlrrs generally shows an 
obvious "bubble structure", consisting of myriad spherical 
low density regions or "holes" in the gas with a large range 
of sizes. 

The origin of these holes has long been attributed 
to stellar feedback by single-age new-born ste llar clusters 
(jWeaver et al.1 Il977l : iMcCrav fe Kafatosl Il987l . and refer- 



ences therein). However, for the LMC it has proven to 
be not at al l evident to correlate HI holes or shells wit h 
Ha emission ( Kim et al ]|l999l : iBook. Chu. fe Gruendlll 20081 ) . 
Holmberg II has similar issues, with Hct not tracing the 
holes, and the stellar ages found therein not corresponding 
well with the kinematical age of the holes ([Stewart et al.l 
I2OOOI : iRhode et all 1 19991: IWeisz et alJl2009h . Studying this 
last galaxy in detail, IWeisz et al.l ([20091 ) propose a multi- 
age model, where HI holes are created by stellar feedback 
from multiple generations of star formation spread out over 
tens to hundreds of Myr. This model is supported by the 
fact that Ha and 24 micron emission, which trace the most 
recent SF, do not correlate well with HI holes, while UV 
emission, which traces SF over roughly the last 100 Myr, 
correlates much better. The concept of a single age for a 
hole is rendered ambiguous. 

In Fig. [11] we show the structure of 3 of our simulated 
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Figure 9. SFHs of our productions runs, see Table [2] The same 
runs as in Fig. [Hare plotted here, comparing spherical and rotat- 
ing models for different masses. Both their SFH (solid lines) and 
evolution of their stellar mass (dashed lines) are shown. 

dwarf galaxies, in a sequence of snapshots taken throughout 
their entire evolution (all shown face on). The projected gas 
density is rendered as the background color (see colorbar), 
and two different age selectio ns of the ste ll ar pop ulation are 
plotted. In accordance with IWeisz et al.1 (|2009h we choose 
these to represent the newest stars (yellow dots, stellar age 
< 20Myr) which would be detected in Ha, and the recent 
stars (red dots, 40Myr < stellar age < lOOMyr) which would 
show up in UV. The gap between the two populations serves 
to provide a clearer distinction between them on the plots. 
These 3 simulations compare a spherical model (205, left 
column), a flattened non-rotating model (265, right column) 
and a rotating model (225, center column), all based on the 
C05 model (see Tables [T] and [2}. The specific snapshot times 
have been selected to represent "interesting" moments in the 
galaxies' SFHs, coinciding with SF peaks or lulls, see Figs. 
[9] and [10] A full, high quality animation, corresponding to 
Fig. [11] can be found online! 



2 HD video: http://www.youtube.com/watch?v=L20WqfMlazo 
YouTube channel of Astronomy department at Ghent University: 
http: / / www.youtube.com/user/AstroUGent 

YouTube playlist with all additional material for this paper: 
http://www.youtube.eom/user/AstroUGent#g/c/EFAA5AAE5 
C5E474D 



Figure 10. SFHs of our productions runs, see Table [2] The runs 
with q = 0.5 from Fig. [8] are plotted here, comparing spherical 
and flattened models for different masses. Similar to Fig. [9] 



4-3.1 Spherical simulations 

In the spherical models gas collapses to the center and forms 
stars that collectively blow out the gas through feedback, 
preferably to one side in a so-called chimney (see snapshots 
at t — 1.17). Over time, the gas cools and re-collapses after 
which star formation can resume again (snapshots at t — 
2.74, t = 5.77 and t = 8.12). This cycle continues throughout 
the entire evolution. There is no significant difference in the 
correlation with the local gas density between the two stellar 
populations shown in Fig.[TT] Both populations are centrally 
concentrated, and so is the gas density. 

Overall there is little small-scale structure: the be- 
haviour of the gas takes place on a large, collective scale. 
This becomes particularly apparent when comparing to the 
rotating model, discussed further on in 14.3.31 An occasional 
small bubble can be spotted in the gas when the galaxy is 
forming stars (e.g. snapshot at t = 4.69, on the upper side 
of the galaxy) . This large-scale behaviour translates into the 
characteristics of the SFHs of the spherical models discussed 
before: large SF peaks separated by quiescent periods. 



4-3.2 Flattened simulations 

The structure of the gas in the flattened dwarf galaxies is 
quite similar to the spherical ones. Large-scale behaviour 
with a centralized structure is still very much the case, which 
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Figure 11. A series of snapshots from the evolution of 3 simulated dwarf galaxies based on the C05 DG model. Left column: basic spherical 
model (205), middle column: rotating model (225), right column: flattened, non-rotating model (265). The bottom of the left 3 columns 
continues on the top of the right 3 colums. Snapshots show rendered gas density (colorbar), new star particles (yellow dots, stellar age 
< 20Myr) and recent star particles (red dots, 40Myr < stellar age < lOOMyr). The label in the top right corner indicates the time in the 
simulation (Gyr), and all galaxies are shown face-on in the x — y plane, except the last 3 snapshots which are edge-on in the x — z plane 
(axes are in kpc). A full, high quality animation can be found online, see text. 
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again can be connected to the discussion and conclusions 
about the SFHs of the flattened galaxies. Small-scale struc- 
ture is not significantly more present than in the spherical 
models, and the previous discussion of the evolution of the 
spherical models is equally valid for the initially flattened, 
non-rotating models. 



4-3.3 Rotating simulations 

The structure of the gas content of the simulated dwarf 
galaxies is noticeably different when adding rotation. There 
is now much more small-scale structure in the gas. A "bubble 
structure" emerges in the gas, caused by the stellar feedback 
of individual star particles (snapshot at t = 0.29) or small 
pockets of star particles (very clear at e.g. t = 6.84). Ap- 
parently, the influence of stellar feedback has become more 
local, and the gas does not exhibit the same global, large 
scale behaviour seen in the rotationless models. 

There is now a very strong difference in the correla- 
tion between the local gas density and the separate stel- 
lar populations. The newest stars are always found in the 
densest regions of the gas, which is not unlogical c onsid- 
ering the star formation criteria ([Valcke et al l l2008h . The 
slightly older stars are much more likely to be found in the 
bubbles or holes because individual groups of star parti- 
cles have had sufficient time to accumulate enough collective 
feedback. T h is all speaks in favor of the mult i- age model of 
IWeisz et al.l (12009) for creating HI holes, and the findings 
of IStewart et al.1 (hooch that young stars (Ra) prefer high 
density HI regions while older stars (FUV) are more likely 
found in low density regions. The idea that UV should be 
a better tracer for HI holes than Ha therefore seems very 
plausible. We can also spot cases of triggered secondary star 
formation, the clearest example being at t — 4.69 where a 
large bubble at the lower left side expands outwards and 
compresses the gas along a rim on the outside of the bub- 
ble, spawning new star formation in this rim. Observational 
evidence for similar events can be found where secondary SF 
is detected in Ha along rims aroun d HI holes ( St ewart et al.l 
l200d : iBook. Chu. fe GruendlH2008T ). 

All this again translates into the SFH characteristics 
we discussed before for rotating dwarf galaxies, where the 
periodicity, or in other words the large scale oscillation, of 
the SFHs from the spherical models is significantly reduced. 
At times when in the spherical and flattened models star 
formation has almost completely ceased, the rotating model 
still shows a significant activity. It continuously forms stars 
throughout the entire simulation. 

As a last point, the SF is also noticeably more spa- 
tially extended than in the spherical and flattened cases. 
Moreover, the spatial extent is quite constant during the 
entire simulation. Stars are always formed throughout prac- 
tically the entire body of the galaxy, while in the spher- 
ical/flattened cases the subsequent SF bursts become in- 
creasingly centrally concentrated. 



4.4 Scaling relations 

Aside from the specific characteristics of individual models 
we discussed above, we also consider the global photomet- 
ric and kinematical scaling relations traced by the simulated 



galaxies. Our main aim is to see how well the general char- 
acteristics of our simulated dwarf galaxies agree with obser- 
vational data of dwarf galaxies as a class. 

In the following, two series of simulations are plotted. 
Firstly, all non-rotating galaxy models, both with spheri- 
cally symmetric and with flattened halos, represented with 
white symbols in Figs. 1121 to 1191 Secondly, all galaxy models 
initially rotating at v ro t = 5 km/s, both with spherically 
symmetric and with flattened halos, are represented with 
grey symbols in these figures. The symbol shapes distin- 
guish the initial halo flattenings: simulations with q = 1 are 
shown with lozenges, q = 0.5 with triangles, q = 0.1 with 
squares. See Tables Hand El 



4-4-1 Half-light radius R e 

The effects of flattening and rotation on the half-light radius 
are clear in Fig. 1121 where R e versus My is shown. Overall, 
at a fixed luminosity, rotation causes R e to increase. The 
initial flattening of the halo does not seem to make a signifi- 
cant difference since both model sequences are quite narrow. 
Only at the high mass end of the non-rotating series does R e 
increase with flattening, the rotating series are unaffected. 
Overall, there does not seem to be a second parameter ef- 
fect. Both series together nicely encompass the observational 
width of the scaling relation. The third series of simulations 
(with a low initial rotation speed of v ro t = 1 km/s) were 
omitted for clarity of the plot. They simply lie between the 
two plotted series, providing us with rotation as a possible 
explanation for the width of the scaling relation. 



4-4-2 Velocity dispersion a 

Fig. [13] shows the stellar central velocity dispersion a versus 
My, projected along the line of sight. We ta ke this to be 
the x- axis, viewing the models edge-on. As in I Valcke et al.l 
(2008), the central velocity dispersion is in general somewhat 
too high. However, the rotating models, having lower veloc- 
ity dispersions, compare favorably to spherically symmetric 
or flattened ones. This decrease of the velocity dispersion in 
the more massive models is tied to the increase of the half- 
light radius in the more massive models. In the rotating 
series the initial flattening also appears to have somewhat 
of an effect, leading to slightly higher velocity dispersions in 
the most flattened cases (squares). This is possibly due to 
the high M/L ratio of these systems, since they have a con- 
siderably higher dynamical mass (for a given stellar mass) 
than their initially less flattened counterparts (see Table [3} . 



44.3 ColorV-I 

The global V — I colors of the models are shown in Fig. 1141 
The rotating galaxies lie a little lower than the non-rotating 
on this plot, meaning these galaxies are slightly bluer. This 
can be understood from their star formation histories (Fig. 
[9}. The strength of the first SF peak is reduced, while at 
later times SF is enhanced with respect to the non-rotating 
case, producing more younger, bluer stars. Otherwise, all 
simulations fall well within the observational range. This 
is, however, not a very stringent test of the models, given 
the fact that the V — I color of an intermediate- age stellar 
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Figure 12. Half-light radius versus V-magnitude. White sym- 
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20041 7wiL04), Ursa Major from 



population is relatively insensitive to metallicity (see next 
paragraph) . 



enrichment are much more diffuse, providing (on average) 
less-met alrich gas for subsequent stellar generations. 



4-4-4 Metallicity 

The metallicity of all dwarf models is too high, especially in 
th e low-mass reg i me. T his problem was already encountered 
bv lValcke et al.l (12008) for the spherically symmetric mod- 
els. Below Mb ~ — 12 mag, the rotating models are more 
metalrich than non-rotating ones whereas above this mag- 
nitude they are less metalrich. An explanation for the low 
mass models can again be found in the respective SFHs (Fig. 
[9])- In the least massive non-rotating models, the large first 
peak in the SFH strongly inhibits further star formation be- 
cause the combined force of the feedback is strong enough to 
severely lower the gas density. Adding rotation reduces this 
first peak and thus also its truncating power, allowing SF to 
proceed continuously and enrich the gas further with subse- 
quent stellar generations. When going to higher masses how- 
ever, the effect and importance of the first peak decreases. 
From the C07 model on, the trend reverses. This is most 
likely due to the strong decrease in central concentration of 
SF and feedback because of rotation (see Figs. ITT1 and l20l at 
the end of this paper), together with the simple fact that 
less stellar mass is produced. SF, metal production and gas 



4-4-5 Surface brightness profiles 

The surface brightness profiles are fitted with a Sersic law 

I(R) = I e y R o) , (7) 

from which the parameters fio and n are plotted and com- 
pared with observational data in Fig. [16] The Sersic index 
n does not differ significantly between rotating and non- 
rotating models. The central surface brightness /xo, on the 
other hand, is consistently lower in the rotating models. This 
is to be expected: with rotation the SF becomes less centrally 
concentrated and more widespread, lowering the central sur- 
face brightness. 



4-4-6 Fundamental plane 

The fundamental p lane ([Bender, Burstein, fe Faberl Il992l : 
iBurstein et al.|[l997h is shown in physical coordinates (Fig. 
fTTj) and in n space (Fig. ITS]). The "vertical" deviation from 
the fundamental plane is shown in Fig. 1191 Except for the 
most massive non-rotating galaxy models (white lozenges), 
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Figure 14. V — / color versus V-magnitude. Symbols as in Fig. 
1121 and the typical error bars for the Fornax Cluster dSph data 
are shown. Observational data as in Fig. [12] 



Figure 15. Metallicity in log 10 (Z [Z©])b (weighed with 
B-ban d luminosity , with = 0.02) versus B-ma gnitude. 

Data :lMateolf l998 (MA98), iNagashima & Yoshiill2004l (NA04), 
iMichielsen et al.ll2007l (MI07). 



which are very compact, most dwarf galaxy models lie sig- 
nificantly above the fundamental plane. These compact non- 
rotating dwarfs have small R e and consequently high mean 
surface brightness within R e (denoted by J e ), making them 
stick out in the side-view of the fundamental plane (Fig. I17p 
and in its m — k,2 projection (Fig. [T8|) . 

Overall, the simulations agree very well with the ob- 
servational trends and luminosity dependent deviations 
from the fundament al plane (observational data taken from 
iBurstein et al.|[l997h . Even the compact non-rotating mod- 
els fall within the observational spread of the relation. 



5 RESULTS / DISCUSSION 
5.1 Evaluation of analysis 

From the previous paragraph, it is clear that rotation has 
a more pronounced influence on the observational proper- 
ties of the simulated dwarf galaxies, quantified by photo- 
metric and kinematical scaling relations, than the flatten- 
ing of the initial conditions. The differences between the 
sequences of rotating and non-rotating models are signifi- 
cantly larger than those between flattened and non-flattened 
galaxies within each sequence. Still, all models fall within 
the range allowed by the data, apart from the problems 
we noted with metallicities being to high. Moreover, despite 
their simplicity, this suite of simulations suggests a possible 
explanation for the widths of the observed scaling relations. 
While mass is the dominant parameter that determines the 
shape and slope of each scaling relation, angular momentum 
could be an important second parameter that determines 
the width of the relations. This will however not be the only 
factor, since external influences such as environment and 
merger history are likely to also have a significant influence 



12 
14 
16 
18 

> 

o 20 
22 
24 
26 
28 

5 

4 

£3 
2 
1 





★ 


DR05 dE 


• LG free dSph 


★ 


GG03 dE/E 


O Perseus dSph/dE 


• 


LG MW dSph 


• Fornax dSph 


o 


LG M31 dSph 






-25 



Figure 16. Sersic parameters versus V-magnitude. Upper panel: 
central surface brightness uo in the V-band, lower panel: Sersic 
index n. Observational data as in Fig. [12] 
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Figure 17. Side-view of the fundamental plane in physical co- 
ordinates. White symbols are our non-rotating galaxies; the grey 
ones are the fastest rotating galaxies, shapes denote initial flat- 
tening as in Fig. fT2l All other points o n the plot are observational 
data, taken from Burstein et al 1ll997l (BU97) and De Rii cke et all 
120051 (DR05). 




Figure 18. The fundamental plane in tt-space. Upper panel: side 
view of fundamental plane (^3, /ci), lower panel : face-on view 
(«2) Observational data from lGeha et al.ll2003l (GE03), other 
symbols as in Fig. [TT] 



here. And we should also note the inherent variance that 
is present in our models, as discussed in section 13.3.11 and 
shown in Fig. [2J 

While the effects of flattening and rotation on the ob- 
served scaling relations are modest, the addition of rota- 
tion has a strong effect on the details of the evolution of 
dwarf galaxies. This is most clearly seen in the properties of 
the stellar populations, e.g. in the metallicity profiles (Figs. 
H and [8]), SFHs fFigs. l9l and [T0|h and overall appearances 
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Figure 19. Deviation from the fundamental plane. Symbols as 
in Fig. [17] 



(Fig. Ill[) . In this respect, rotating models are qualitatively 
quite distinct from non-rotating ones, independent of ini- 
tial flattening: rotating models have continuous SFHs with 
wide-spread SF while non-rotating models have "breathing" 
SFHs with centrally concentrated SF. Observationally, this 
leads to rotating models having flat metallicity profiles while 
non-rotating models show pronounced negative metallicity 
gradients. 



5.2 Mechanism 

Within the sequence of non-rotating models, flatttening of 
the initial conditions appears to have very little effect on 
the models' properties, which is of importance especially 
when considering the metallicity profiles. It therefore seems 
doubtful that the fountain mechanism is very important 
for dwarfs. Still, if there are large feedback driven out- 
bursts of gas, they tend to be aligned preferentially along 
the minor axis. But the expelled enriched gas does not fall 
back onto the galaxy. This is most likely because of the 
shallow potential wells of dwarf galaxies, and because the 
remaining cold gas is simply "in the way". Another im- 
portant argument against the fountain mechanism is the 
actual flattening of dwarf galaxies, both in observations 
and in our simulations, as we discussed in section 13.3.31 
Dwarf g alaxies simply are not likely to occur with very fla t 
shapes (|Sanchez-Janssen, Mendez-Abreu, &; Aguerrll2010h . 
Their flattenings are not comparable to those of massive 
spiral galaxies (q ~ 0.2), they are on average much thicker 
((q) ~ 0.6). This often makes it difficult to even speak of a 
"disk" in the context of dwarf galaxies. Therefore, while the 
fountain mechanism might be very relevant in the domain 
of large spiral galaxies, with much deeper potential wells 
and much flatter shapes, it does not appear an important 
mechanism in dwarf galaxies. 

Rotation, on the other hand, leads to important quali- 
tative and quantitative changes in the SFHs of dwarfs. The 



16 J. Schroyen et al 




time [Gyr] time [Gyr] 



Figure 20. Star formation density in M0/Gyr/kpc 3 for spherical, non-rotating (left column) and rotating models (right column). The 
total mass of the models increases from the top down as indicated in the figure. SF density is plotted in color-code versus time on the x 
axis and radius on the y axis. 



consequences of the addition of angular momentum are the 
following: 

(i) Gas will spiral inward, instead of falling straight to 
the center. There is a " centrifugal barrier" preventing the 
gas from collapsing to a dense central region. 

(ii) Since the gas density is much more smeared out, so is 
the star formation. The density criterion (see section [2.1. 1[) 
for star formation is now reached in a much larger region of 
the gas, so that star formation will occur throughout prac- 
tically the entire body of the galaxy. This is evident in Figs. 
1111 and [201 star formation is consistently more spatially ex- 
tended in comparison with non-rotating models. 

(iii) This naturally produces more spatially homogeneous 
stellar populations. Therefore the gas is enriched much more 
homogeneously across the entire galaxy, explaining the flat 
metallicity profiles in Fig. [7] 

(iv) Where there is star formation, unavoidingly there 
will also be stellar feedback. Since the former is smeared 
out across almost the entire galaxy, so is the latter. The 
supernova feedback now being less centrally concentrated, 
this leads to much less pronounced large-scale collective be- 
haviour of the gas. The effects of feedback are now more 
local. This has two distinct but related effects: 

• The supernovae combine their energy locally on a 
smaller scale, and produce low-density holes in the gas, 
instead of collectively blowing out the gas and lowering 
the global gas density after a large centralized star forma- 



tion event. This hole or bubble structure is clearly visible 
in Fig. Illland is discussed in section l4~3l 

• This can also be linked to our findings concerning 
the SFHs in section 14.21 and Fig. [9] Since the gas does 
not collectively blow out due to feedback, star formation 
will not shut down completely across the entire galaxy, 
because only locally the density criterion for star forma- 
tion is not satisfied (in the feedback holes, Fig. 1 1 1 1) . This 
is also seen in Fig. [20] Collective behaviour - i.e. large- 
scale oscillations in the SFR - is diminished, leading to 
more continuous, less variable SFH. The "breathing" SF, 
typical of non-rotating models, is largely absent. 

The density criterion mentioned in section [2] and here 
in point (ii) is an important element of our models. We 
should note that we em ploy a treshold of 0.1 cm -3 , while 
iGovernato et al.1 (|2010h suggest the usage of a treshold of 
100 cm -3 , reflecting more realistically the conditions of real 
star-forming gas clumps. We do not expect this to qualita- 
tively change the proposed mechanism however: the higher 
treshold will take the gas longer to reach it when collaps- 
i ng, but the exte nsion of the cooling curves below 10 4 K 
(Ma io et al.ll2007h will cause the gas to collapse easier and 
on smaller scales. These effects might not cancel each other 
out, but the relative influence of added rotation, as dis- 
cussed above, will remain qualitatively similar. Perhaps on 
small scales the chemical homogeneity will be less, but on 
large scales rotating galaxies will still be chemically homo- 
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geneous. The inclusion of the high-density treshold and the 
extra cooling will be the subject of further research. 



5.3 Galaxy types 

As already mentioned in paragraph 14.41 our dwarf galaxy 
models agree quite well with the observed the scaling rela- 
tions of early- type galaxies. However, since the model galax- 
ies still contain gas and have ongoing star formation at the 
end of the simulation, they should be classified as late- type 
dwarfs. 



5.3.1 dlrrs? 

The non-rotating and slowly rotating models, both flattened 
and non- flattened, are characterized by 

• centrally concentrated gas distribution; high central 
density 

• low specific angular momentum 

• strong stellar population gradients 

• bursty or episodic SF 

• centrally concentrated SF 

• large-scale feedback driven outflows and a largely fea- 
tureless ISM. 

The fast rotating models, both flattened and non- flattened, 
are characterized by 

• spatially extended gas distribution; low central density 

• high specific angular momentum 

• small stellar population gradients, if any 

• continuous SF 

• small star forming regions, scattered across the galaxy 

• turbulent ISM with distinct feedback driven holes. 

dlrrs are known to have a more extensive and less 
centrally concentrated gas distribution that other gas-rich 
dwarf galaxy types (e.g BCDs), and also a relatively high 
specific angular momentum. Chemical homogeneity is a 
general trait of dlrrs, both in th e ir gas and their stellar 
content (jTolstov. Hill fe Tosil [20091 : iKobulnickv fe Skillmanl 
19971). From the r e view of dwarf galaxy properties in 
Tolstoy. Hill fe Tosil (|2009h and the extensive work of 
Dolphin et al.l (|2Q05h , using CMP analysi s to reconstruct 



dwarf galaxy SFHs jAparicio et all |l996l: iDolphinl Il997l 
l2002l : lTolstov fe Sa ha 1996: lTosi et al.lll99lh . it is clear that 



dlrrs generally have a "continuous" SFH without quiescent 
periods without SF. The characteristic gas structures of 
dlrrs have already been discussed in section 14.31 

From this short overview of the observed properties of 
late- type dwarfs, it is clear that our fast rotating models 
resemble dlrrs, at least qualitatively. But our non-rotating 
and slowly rotating models do not, although they do also 
still contain gas and show ongoing (periodic) star forma- 
tion. Angular momentum, it seems, invokes different star 
formation modes in dwarf galaxies. It differentiates between 
centralized/bursty and extended/continuous star formation, 
and all dwarf galaxy properties connected with this which 
are mentioned above. Although quantitatively not compa- 
rable to our models, it is worth mentioning BCDs and 
their differences with dlrrs. They too are gas-rich late- type 
dwarf galaxies, but have a lower specific angular momentum 



and much more con centrated gas distribution (|van Zee et al.l 
l200ll : Ivan Zed 20021). BCDs also show substantial color gra- 
dients ([van Zed l2002h , indicating chemical inhomogeneity, 
and by definition have bursting SFHs. 



5.3.2 Conversion of late-type dwarfs to early types 

Since internal processes such as supernova feedback are 
not capable of removing the gas from a dwarf galaxy, 
we turn to external or environmental processes, e.g. tidal 



g. tida l 
l2006h . 



stripping and ram pressure stripping ([Maver et al. 
Ram-pressure stripping is able to remove a large frac- 
tion of the gas and leaves the structure and kinemat- 
ics of the stars relatively undisturbed, thus preserving 
any pre-existing stellar population gradients and rotatio n 
(jGrebel et al.ll2003l : iMarcolini. Brighenti. fc D'ErcoldliooiSh . 
Tidal interactions can cause violent reactions in dynamically 
cold t hin-disk dwarf galax ies and can significantly disturb 
them ([Maver et al ] |2001al lbh. However, the majority of the 
dwarf late- type population is quite round, with mean axis 
ratio (q) ~ 0.6. In such gala xies, tidal interactions wreak 
much less havoc (|Valckdl2010h . 

We therefore argue that it is possible to convert late- 
type dwarfs into early- type ones inside a cluster environment 
by removing their gas and halting SF without significantly 
altering their structural and kinematical properties. So the 
rotation which is present in dlrrs can be preserved in their 
dE descendants along with the stellar characteristics con- 
nected with rotation (metallicity profiles). 



6 CONCLUSION 

The centrifugal barrier mechanism formulated in section [5] is 
able to combine all our findings we discussed in the analysis 
into one coherent picture, emphasizing the importance of 
rotation in dwarf galaxy behaviour. 



6.1 Metallicity profiles 

Our int erest in this subject was triggered initially by the 
finds of IKoleva etHI <l2009h , who found that dwarf early- 
type galaxies without stellar population gradients were also 
the fastest rotating ones. We conclude from our simulations 
that (in isolation) rotation, or the absence thereof, is indeed 
a key factor in creating stellar population gradients. The 
"fountain mechanism" does not seem relevant on the scale 
of dwarf galaxies, and our simulations clearly indicate that 
the geometry or flattening of a dwarf galaxy does not have 
any significant influence: pressure-supported, non-rotating 
systems behave very much alike, independent of flattening. 

We therefore propose the alternative "centrifugal bar- 
rier mechanism" in section 15.21 which explains the existence 
of flat metallicity profiles as a natural consequence of its 
rotation. 



6.2 Angular momentum as second parameter 

We suggest angular momentum as being a crucial second 
parameter in determining the appearance and evolution of 
dwarf galaxies, with the total galaxy mass being the promi- 
nent first parameter. While our simulations are admittedly 
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very idealized and cannot purport to paint a cosmologically 
up-to-date picture of dwarf galaxy formation, they have the 
enormous benefit of allowing us to unambiguously identify 
the influence of individual parameters, such as angular mo- 
mentum. 

We have shown that rotation has a significant impact on 
the stellar populations of dwarf galaxies. And in the same 
vein we can say the opposite for dwarf galaxy flattening, 
which shows no significant influence in our simulations, and 
thus is less likely to be a major player in dwarf galaxy evo- 
lution. 

6.3 Making dlrrs 

We find that without rotation, it does not seem possible 
to qualitatively produce "typical" dlrrs with spatially ex- 
tended SF, continuous SFHs, a turbulent ISM with low- 
density holes, and most importantly with chemical homo- 
geneity throughout its body of gas and stars. Non-rotating 
models do not display any of these characteristics (having 
centralized SF, bursty SFHs, featureless ISMs and metallic- 
ity gradients). Angular momentum appears to differentiate 
between bursty and continuous star formation modes. 
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